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Abstract. The nodal lines of random wave functions are investigated. We 
demonstrate numerically that they are well approximated by the so-called SLEg curves 
which describe the continuum limit of the percolation cluster boundaries. This result 
gives an additional support to the recent conjecture that the nodal domains of random 
(and chaotic) wave functions in the semi classical limit are adequately described by 
the critical percolation theory. It is also shown that using the dipolar variant of SLE 
reduces significantly finite size effects. 
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1. Introduction 

In 1977 Berry conjectured pQ that the wave functions of chaotic quantum systems can 
statistically be described by Gaussian random functions with a spectrum computed from 
the ergo die average. 

For example, the wave functions for two-dimensional billiards without magnetic 
field ty(x) have to obey the equation 



and, say, the Dirichlet boundary conditions. 

From general considerations it follows that such functions can be represented as a 
formal sum over elementary solutions of e.g. 



where J m (x) are usual Bessel functions and C_ m = C^. 

The Berry conjecture signifies that for chaotic billiards like the stadium billiard, real 
and imaginary parts of coefficients C m with non-negative m are independent identically 
distributed Gaussian random variables with zero mean and variance computed from the 
normalization. 

The beauty of this profound conjecture lies first of all in its simplicity and generality. 
On the other hand, it represents chaotic wave functions as completely structureless 
and 'uninteresting' objects which to some extent was the reason of relatively few 
investigations of the chaotic wave functions. 

In 2002 Smilansky et al. j2] stimulated a renewal of interest in this problem. Instead 
of considering the coefficients of an expansion as in (j2J) these authors looked for the 
nodal domains of wave functions, i.e. the regions where a function has a definite sign. 
They have also observed that for chaotic wave functions the nodal domains have a rich 
unexpected structure worth investigating in details. 

In [H] on physical grounds it was conjectured that the nodal domains of random wave 
functions (and, consequently, of chaotic wave functions) can be adequately described 
by a critical percolation model (see e.g. pj). Due to the universality of the critical 
percolation it is unessential what special percolation model is considered. All of them 
lead to the same critical exponents as well as to other universal quantities. In jH] it 
was checked numerically that many predictions of the critical percolation model, such 
as the number of connected domains, their area distribution, cluster fractal dimensions 
etc. agree very well with corresponding quantities computed from the nodal domains of 
random (and chaotic) wave functions. In [3] it was demonstrated that the level domains 
(the regions where a function is bigger than a certain non-zero value) are well described 
by a non-critical percolation model. These and other investigations strongly suggest that 
the critical percolation model is applicable for the description of the nodal domains of 
random and chaotic wave functions. Though it sounds physically quite natural and can 
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be confirmed by a careful application of Harris' criterion a rigorous mathematical 
proof was not yet found and further numerical verifications are desirable. 

The purpose of this paper is to check another prediction of the percolation model 
namely that the boundaries of the percolation clusters are described by what is called 
SLEg curves (see below). This statement was proved in 2001 for the critical percolation 
on the triangular lattice by Smirnov [7^, and it is widely accepted that it remains true 
for all critical percolation models. For the nodal domains of random functions it leads 
to the conclusion that the nodal lines (curves where a function is zero) have to be also 
described by SLE 6 curves. 

The plan of the paper is the following. In section El for completeness we give an 
informal introduction to the SLE curves. This name is attributed to two-dimensional 
self-avoiding curves generated by an one-dimensional Brownian motion with zero mean 
and the variance linear in time with a coefficient of proportionality equal to a real positive 
number k. In 2000 Schramm proved that if the different parts of a random self-avoiding 
curve are independent and the curve itself is in a certain sense conformally invariant then 
it belongs to SLE curves with a certain value of k. Smirnov's result [7] means that the 
percolation boundaries are generated by the Brownian motion with k = 6. In sectional 
we present the results of direct numerical calculations of corresponding Brownian-like 
generating curves for the nodal lines of random wave functions and found that they 
are well described by SLEg. The inevitable drawback of numerical calculations is that 
one always deals with curves of finite size. But the theorem about the relation between 
the percolation boundaries and SLEg is valid only for infinite curves. To decrease such 
a finite size effect we use in section 0] a different version of SLE called dipolar SLE 
[TO] . In this approach one first conformally transforms a given region into, say, an 
infinite band and then uses a suitably modified SLE equation. By this method we 
numerically demonstrate that the nodal lines of random wave functions agree well with 
the SLEg description thus once more confirming the conjectural relation between the 
critical percolation and the nodal domains of random wave functions. In the appendix 
the numerical algorithms used in calculations are shortly discussed. 

To a certain degree our work was stimulated by [5] where it was numerically checked 
that the nodal lines of the vorticity field in two-dimensional turbulence are close to SLEg 
curves. In the recent publication ^JJ Keating, Marklof and Williams have demonstrated 
that the nodal lines for a perturbed cat map are well described by SLEg. They also 
performed [lSj numerical calculations for the nodal lines of random wave functions which 
lead to the same conclusions. Our results are in complete agreement with their findings. 

2. Schramm-Loewner Evolution 

Two-dimensional self-avoiding curves appear naturally in many physically important 
problems. But analytically imposing the condition of self-avoiding is not simple. In 1923 
Loewner proposed [S] to describe such curves by conformal transformations which map 
a domain with a simple curve growing from the boundary to another simply connected 
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domain without the curve. In the simplest setting one considers the upper half plane H 
with a simple curve C and looks for a conformal map which transforms the upper half 
plane minus the curve, M\C, to H itself (see figure HJ). 
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Figure 1. Loewner's evolution. The upper half plane of the variable z with a cutoff 
along a simple curve is transformed by a conformal map w — g(z) to the whole upper 
half plane of the variable w. 

The main point of this approach is the well known Riemann theorem according 
to which any simply connected region can be conformally mapped to another simply 
connected region and, inversely, if one applies a conformal transformation to a simply 
connected region the result will be also a simply connected region. 

Such conformal transformations are not unique and depend on 3 free parameters. 
For example, the upper half plane EI is transformed to itself by 3 parameter group of 
fractional transformations 

az + b 

9{Z) = —d (3) 
with real a, b, c, d. 

To fix uniquely the conformal map g(z) it is convenient to impose the so-called 
hydrodynamic normalization by fixing the behaviour of the map at infinity: 

2t 

g(z) = z H h • • • when \z\ — * oo . (4) 

z 

Parameter i is a characteristic of the initial region called capacity. To stress its 
importance we shall denote the map g(z) with asymptotic (0} by gt(z) showing explicitly 
the dependence of t. This quantity can be also called 'conformal time' as it fulfils the 
additivity property which follows directly from (J3| 

9t(9s) = 9t+s • (5) 

The simplest example of a curve is a vertical slit in the upper half plane. If (£, h) are 
coordinates of the top of the slit (see figure EJ) the transformation 

w(z) = Z+y/(z-tf + 4t (6) 

with t = h 2 /4 maps the upper half plane of z with the slit to the upper half plane of w 
without the slit. 
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Figure 2. Transformation of a vertical slit to the upper half plane by map 0. 



Interpreting t as the time permits to derive the equation for gt(z). From © it 
follows that 

9t+st(z) = gst (9t(z)) ■ (7) 

For small St any line can approximately be considered as a vertical slit growing from 
the image of the tip of the curve under g t which we denote by £ t (see appendix). Now 
(0) states that 

gst{z) — >z + -. (8) 

z — qt 

Combining two last equations one concludes that gt(z) as function of t obeys the 
following differential equation called Loewner's equation 

dgt(z) _ 2 
dt g t (z) - & 
with initial condition 



(9) 



go(z) = z . (10) 

The function £ t is called the 'driving' or 'forcing' function. The curve itself is called the 
'trace' and can be computed from the relation 

z(t)=gr 1 &). (11) 

The advantage of Loewner's equation @ and (jllj) is that for any smooth ( different iable) 
function £ t it will produce a two-dimensional self-avoiding curve and, inversely, any self- 
avoiding curve will give a smooth function £ t . 

The next important step was done by Schramm. In 2000 he proved that random 
self-avoiding curves with the following properties: 

• conformal invariance, 

• statistical independence of different parts, 

• reflection symmetry 
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are generated by the usual one-dimensional Brownian motion with zero mean and the 
variance linear in t 

& = y/nB t , < B t >= , < B t B s >= \t-s\. (12) 

Such curves only depend on one real positive parameter k and are called stochastic (or 
Schramm) - Loewner evolution curves or shortly SLE K . 

As Brownian curves are not different iable, SLE K curves may have self-touching 
points. Depending on the value of k they are divided into three phases (see e.g. [T5]). 
For < k < 4 the traces are simple curves, for 4 < k < 8 they can have double points 
or hit the real axis, and for k > 8 they fill the entire domain. 

SLE is a very powerful tool to study rigorously the scaling (continuum) limit of 
different discrete models (see e.g. jTHI and |14j ) . The most important for us is the 
Smirnov's result [7j: the boundaries of clusters in the critical percolation (on a triangular 
lattice) converge in the continuum limit to the traces of SLEg. As the nodal domains 
of random wave functions are conjectured to be described by critical percolation [3] 
it means that the nodal lines of random wave functions in this limit have to be also 
described by SLE 6 curves. The verification of this statement is the main purpose of this 
work. 

3. The nodal lines of random wave functions 

The nodal lines are lines where a real function of two variables is zero, ty(x,y) = 0. 
They can, in principle, be computed numerically for any function ty(x,y). For random 
wave functions (J2J) the nodal lines are quite complicated. An example of a nodal line 
is presented in figure 01 To demonstrate the complexity and self-avoiding character of 
such a curve, its small portion enclosed by a rectangle is enlarged in the inset of this 
figure. 
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Figure 3. A nodal line of a random wave function. The dashed line indicates the 
absorbing half-circle. Insert: Magnification of the indicated rectangular region. 
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As all pictures of nodal lines are necessarily finite, one has to decide what to do 
with a nodal line which hits the boundary. There are two main types of boundaries: 
reflecting and absorbing. If a growing line touches a reflecting boundary it follows the 
boundary in the increasing direction till the next appearance of the sign change and 
then continues along a new nodal line. This can be achieved automatically by imposing 
the sign function outside the boundary. If a boundary is absorbing, one simply stops 
any line which hits it. In figure El the horizontal line is a reflecting boundary but the 
half-circle indicated by the dashed line is an absorbing one. 

In general the properties of SLE K curves inside a region do depend on chosen 
boundary conditions. But for SLE 6 lines it is proved (see e.g. [IB]) that they feel 
the boundaries only when they hit them so the boundary conditions are unessential 
for the investigation of the local properties of the percolation cluster interfaces and, as 
conjectured, of the nodal lines of random wave functions. 

2 | ■ , ■ , 1 
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Figure 4. Forcing function corresponding to the nodal line of figure 01 

Having generated a line as in figure 01 the next step is to find the conformal map 
which transforms it to the upper half plane. There are many good algorithms which 
permit to do it numerically (see ^3] and references therein). We used the so-called 
geodesic algorithm 15^ in which a small segment of a line to be mapped is approximated 
by an arc of a geodesic circle perpendicular to the real axis and passing through the tip 
of the segment (see appendix). This algorithm is quite stable and is easy to implement. 

In figure 0] we plot the forcing function which corresponds to the nodal line of 
figure H3 and in figure El fine details of this function are given. 

We have numerically computed 2248 different realizations of the random wave 
function (J2J) with k = y/E = 100. For each realization we have found the nodal lines 
which were stopped when they hit the half-circle as in figure El 
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Figure 5. (a) Magnification of the small rectangular region in figure 0] (b) 
Magnification of the rectangular region in (a). 



Then corresponding forcing functions, £t(j), were calculated using the geodesic 
algorithm. To check that these functions are well described by the Brownian curves we 
calculated mean value £j and variance cx 2 (t) for each t from the usual relations 

N N 

& = af £ ' °*® = n £ (6(j) " & a • (13) 

In figure IH1 the dependence o~ 2 (t) is plotted for all t. It is clearly seen that it is not linear. 
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Figure 6. Variance of forcing functions for random waves inside the half-circle. 
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This is not surprising as the linear behaviour is predicted only for infinite curves. 
Curves closed to boundaries which do not exist in SLE require a special treatment which 
we will discuss in the next section. For short time when curves do not feel the boundary 
the variance is close to be linear. 

In figure U\ the variance and the mean value are presented in a short interval of t 
where cr 2 (t) is approximately linear. The solid line in this figure is the best quadratic 
fit to the variance 

y = -0.003 + 6.05t- lO.Ot 2 . (14) 




0.00 0.02 0.04 

t 



Figure 7. Variance (circles) and mean value (dashed line) in a short range. Solid line 
is the best quadratic fit (|14ll . 

The most important for us is the value of the linear term, 6.05. It is quite close to the 
pure percolation value k = 6. The mean value of the forcing functions (~ 0) also agrees 
with the SLE value. 

Another important prediction of the SLE 6 description of the percolation model is 
that the distribution of values of £t with fixed t has to be Gaussian with variance 6t. In 
figure |H1 this distribution is presented for different values of t. When the abscissa axis 
is rescaled by v6i all data are quite close to the Gaussian curve with unit variance. 
Similar results were obtained when instead of the half-circle we consider nodal lines 
inside a rectangle. 

The agreement of our results with SLE 6 predictions confirms that the nodal lines of 
random wave functions are well described by the percolation model, as was conjectured 
in @j. 
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Figure 8. Probability law for random wave forcing functions. The abscissa axis is, 
as it follows for SLEg, u = £ t /y/6t. The solid line is the standard Gaussian with unit 
variance: P(u) = exp(— u 2 /2)/\/2tt. The points indicated by • , A, and □ correspond 
to, respectively, t = 0.02, 0.04, and 0.06. 



4. Dipolar SLE 

Though the results of the previous section are in a good agreement with SLEg, let us 
notice two facts. First, the approximately linear increase of the variance with time 
is observed only for a very small time interval. Second, even within this interval a 
quadratic term in t always exists and is quite large (cf. (|14|)). Both of these drawbacks 
are related with inevitable finite size effects and can be reduced by the using of another 
variant of SLE. 

The situation is analogous to the Brownian motion in a finite interval. It is well 
known (see e.g. ^Bj) that the probability for an one- dimensional Brownian particle 
starting from xq to arrive to the point x at time t in the whole space is given by the 
Green function of the diffusion equation 

G ^-v^M-^f) (15) 

where D is a diffusion coefficient. Standard relations < x >= xq and < (x — xq) 2 >= Dt 
are simple consequences of this expression. 
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For the Brownian motion inside a finite interval such probability depends on the 
boundary conditions. There are two main types of boundary conditions: reflecting and 
absorbing. When a particle hits a reflecting boundary it reflects back, when it touches 
an absorbing boundary it stops. It is known (see e.g. |16j ) within the formalism of 
the diffusion equation that a reflecting boundary gives rise to the Neumann boundary 
condition and an absorbing one corresponds to the Dirichlet condition. When both 
boundaries of an interval of length h are of the same type one gets that the Green 
function of this interval is 

oo 

G±(x,x ;t)= ^ (Go(x + 2mh,x ;t)±G (-x + (2m + l)h,x ;t)) (16) 

m=— oo 

Here + (resp. — ) denote Neumann (resp. Dirichlet) boundary conditions imposed at 
points ±h/2. 

In figure El the mean variance 

fh/2 

a 2 {t)= I (x-x) 2 G±(x,x ;t)dx (17) 

J-h/2 

is plotted as a function of t for the case xq = and D — 1. It is clearly seen that 
the long time behaviour depends on the boundary conditions and that the variance is 
linear in t only for very small t . These curves are similar to figure EH where the linear 

0.4 
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Figure 9. Variance for the one-dimensional Brownian motion with absorbing 
boundaries (solid line) and reflecting ones (dashed line). 

behaviour was observed only for t < 0.05 (cf. figure EJ). We estimated that in the best 
cases this interval corresponds to, roughly speaking, only 1/4 of the total curve length. 

To take into account larger parts of the curves and to check the Brownian-type 
behaviour for a longer time interval it is required to construct formulae similar to (|16j) 
for SLE. 
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Figure 10. (a) Rectangular lattice inside the half-circle, (b) Its image under the map 
(H3- 



The simplest case corresponds to the reflecting boundary conditions. In this setting 
one considers a region and random curves (nodal lines) which start at one fixed point, 
Zo, of the boundary and end at another fixed boundary point, z\. For the nodal lines 
this can be achieved e.g. by imposing that one boundary arc from zq to z\ is positive 
and the second one is negative. This type of processes is called chordal SLE (see e.g. 
[13]). By definition, it can be reduced to the standard SLE from to infinity (cf. 0) by 
a conformal transformation which transforms the given region to the upper half plane 
in such a way that the point Zq is mapped to the origin and the point z\ to infinity. 

Here, we would like to use a different variant of SLE which corresponds to a region 
with two boundary arcs restricted by points z + and z_. One arc is assumed to be 
a reflecting boundary and the other is an absorbing boundary The random curves 
emerge from the point z on the reflecting boundary and are stopped when they hit 
the absorbing arc. By a conformal transformation our region can be transformed to the 
standard strip § 

§ = {z G C,0 <Im z < vr} (18) 

in such a way that points z_, z + , and zq are mapped to — oo, +oo, and respectively. 

In [TO] it was shown that the SLE re process which joins to a point of the line 
Xm z = Ti is described by the following Loewner-type equation 

M^l = ? 

dt tanh(gf t (z) - £t) 

where go(z) = z and the forcing function, £ t , is as above the standard Brownian motion 
(112)1 with variance nt. 

To construct the dipolar SLE numerically from the data used in the previous section 
we first transform the nodal lines inside the half-circle, as in figure El to the standard 
strip (|18j) by the following transformation 

F(z) = ln[(L + z) 2 /(L - z) 2 } (20) 

where L is the circle radius. 

To visualize this mapping we present the images of a rectangular lattice inside the 
half-circle in figure HOI Note the strong deformations of the regions close to the absorbing 
boundary. 
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The first step of calculations consists in the transformation of every nodal lines 
inside the half-circle to lines inside the chosen strip (|1B|). Such an example is illustrated 
in figure HU 

4 | . 1 . 1 . 1 . 1 . 1 




-5 -2.5 2.5 5 

Figure 11. Image of the nodal line of figure under the map l|20() . The dashed line 
is the absorbing boundary. 

Then it is necessary to find the dipolar conformal transformation which map the 
strip with a line to the strip itself (cf. In the appendix a simple algorithm for such 

a mapping is briefly discussed. Using it we compute numerically the forcing function, 
£t for each nodal line and calculate its statistical properties. In figure EH the variance 
(TT3J) is plotted for all curves and in figure H31 the region of linear increase of cr 2 (t) is 
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Figure 12. Variance of forcing functions for the nodal lines inside the strip. 
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magnified. The solid line in this figure indicates the best quadratic fit 

y = 5.92t-0.103t 2 . (21) 

The coefficient 5.92 is to compare with the percolation theory value of k = 6. Other 
statistical characteristics are also in a good agreement with the SLE 6 predictions. 

Contrary to (fT4*|) the quadratic term in (|21|) is small. It is mainly related with 
discretization errors and it decreases when more points along a curve are taken into 
account. The non-linear behaviour of the variance in figure H~2l for large t is connected 
with large errors due to the stretching of lines close to the boundary as it is evident from 
figures El and ^2 We checked that the interval t < 8 in figure EH where the variance is 
linear in t corresponds approximately to 92% of the total curve length. It means that 
the dipolar SLE 6 is a good description of the nodal lines practically for the entire curves. 
Exceptions are the small line parts close to the absorbing boundary which require better 
appr oximat ions . 




2 4 6 8 

t 

Figure 13. Variance (circles) and mean value (dashed line) for 2252 realizations inside 
the strip. Solid line: the best fit i|21[l. 

We have also investigated nodal lines inside a rectangle [—L/2; L/2] x [0;/] which 
can be conformally mapped into the strip (JT5J) by the following transformation 

F{z) = - ln[(p(z + L/2) - p(L))/(p(L/2) - p{L))\ (22) 

where p(z) is the Weierstrass elliptic function with periods 2L and 2il (see e.g. |T7|). 

The results in this case are similar to the ones above and all agree well with the 
percolation model. 
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5. Conclusions 

In summary, we investigated the nodal lines of two-dimensional random wave functions 
by numerically computing one-dimensional forcing functions of Loewner's evolution. We 
demonstrate that the later are well described by a Brownian motion with zero mean 
and variance growing linearly in time. The coefficient of proportionality is close to the 
value 6 predicted by the percolation model. Our results give an additional support to 
the conjecture that the nodal domains of random wave functions in the scaling limit 
are described by the critical percolation. We show that using dipolar SLE reduces 
significantly finite size effects. 
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Appendix 

Consider a simple curve which starts on the real axis and grows inside a region of the 
upper half plane. The goal of the numerical algorithms below is to find a forcing function 
in (jnj) corresponding to a conformal map which transforms the region cut along the curve 
into the region itself. 

Let (x , 0) be coordinates of the first point of the curve which we assume belongs to 
the abscissa axis and let (xi, yi) be coordinates of the next point closest to the first. The 
simplest numerical method consists of assuming that the forcing function is a piecewise 
constant function. The map with a constant forcing function 5£ for < t < 5t is given 
by © as 



The parameters <5£ and 5t are obtained from the condition that the top of the slit 
coincides with the point 5z = x\ — xq + iyi, which yields 



Then we transform all points of the curve except the first one by (|23|) with these values of 
parameters and renumber the resulting points. One repeats this process till the whole 
curve will be transformed. In jTHj it is proved that such an algorithm converges for 
sufficiently small 5tj and Sxj. 

In section we used the more refined geodesic algorithm ^Hj in which one 
approximates a small part of the curve between two points, and Sz, by the geodesic 




(23) 




(24) 
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arc (a circle perpendicular to the real axis) passing through these points. Normalizing 
the map given in ^5] such as to obey the convention (jjj) one gets the expression 

& 2b 3 + 36c 2 

~ VVT^^ibz/ib - z)) 2 + c 2 -b 2 -c 2 + 2(6 2 + c 2 ) 

where 

b _ \Sz\ 2 ^ c _ \Sz\ 2 _ (26) 
TZeSz ' ImSz 

Direct calculations give that the time corresponding to 5z and the value of the forcing 
function in this point are 

5t = - (1m Szf + - {Tie Szf (27) 
4 8 

and 

6£ t = lne6z. (28) 

When a point moves along a geodesic circle of diameter b the forcing function changes 
as follows 

6= ^ (29) 

and is practically linear for short time. Eqs. ()25 |) -()28 |) define the geodesic algorithm. 

To do some computations following dipolar SLE in Section |U we use a piecewise 
constant approximation for the forcing function in ()19j) . In this case one has the exact 
solution [10J 



cosh (~{g t {z) - SO) = e 5t/2 cosh (l(z - 6£)\ . 



(30) 



Finding St and S£ from the condition that the curve tip coincides with Sz one obtains 

S£ = TZeSz , exp ^— ^St^j = cos ^-Imfcj . (31) 

These expressions are the dipolar analog of Eqs. (|2"4"j) and they permit to construct the 
simplest algorithm of numerical calculations for the dipolar case. 
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